Detecting level crossings without looking at the spectrum 
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In many physical systems it is important to be aware of the crossings and avoided crossings which 
occur when eigenvalues of a physical observable are varied using an external parameter. We have 
discovered a powerful algebraic method of finding such crossings via a mapping to the problem 
of locating the roots of a polynomial in that parameter. We demonstrate our method on atoms 
and molecules in a magnetic field, where it has implications in the search for Feshbach resonances. 
In the atomic case our method allows us to point out a new class of invariants of the Breit-Rabi 
Hamiltonian of magnetic resonance. In the case of molecules, it enables us to find curve crossings 
with practically no knowledge of the corresponding Born-Oppenheimer potentials. 

PACS numbers: 03.65.-w, 31.50.Gh, 31.15.-p, 32.60.+i 
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Curve crossing is a universal phenomenon with exam- 
ples in many branches of both pure and appUed sciences. 
It is responsible for electron transfer in proteins hi, un- 
derlies stability analysis in mechanical engineering |3 and 
determines Nash equilibria |3- Mathematically, it arises 
naturally in the subject of algebraic geometry |j|. In 
physics, the discovery of a crossing between two energy 
levels, for example, signals degeneracy in the energy spec- 
trum and the existence of an underlying symmetry in the 
problem. Since the initial pioneering works by Hamil- 
ton [5j, von Neumann and Wigner [6j, and Landau [TJ 
and Zener |3] , many interesting physical phenomena have 
been associated with crossings. These include Berry's 
phase in adiabatic quantum mechanics |9j, the hidden 
symmetries of the Hubbard model in condensed matter 
physics llOl. a nd the onset of quantum chaos in nonlinear 
dynamics |lll |. In contemporary atomic and molecular 
physics an avoided crossing that has come into intense 
focus is the Feshbach resonance. Such a resonance pro- 
vides unique experimental control over interactions in a 
quantum deg enerate atomic gas, realizing collapsing con- 
densates [13, ultracold molecules [l3j and the crossover 
from BEC to BCS pairing in degenerate Fermi gases [1J| . 

In this Letter we present a versatile theoretical tech- 
nique for detecting the presence of level crossings in phys- 
ical systems. Our technique is algebraic in nature, as the 
underlying Hamiltonian (or other physical quantity) can 
often be represented as a parameter-dependent matrix. 
Our method is powerful and general, as wc extract essen- 
tial information simply from the matrix elements with- 
out having to compute the spectrum, and remarkably, 
works even in cases where the Hamiltonian is not fully 
known. We demonstrate our technique on atoms and 
molecules in a magnetic field. Apart from being interest- 
ing in themselves, these physical systems are of relevance 
to Feshbach resonances, a topic of current experimental 
and theoretical concern [l3|. We present specific results 
for the atom pair ^^Na-*^Rb, which is of current inter- 
est as a candidate for Feshbach resonances [l5|, but our 
analysis can be readily extended to other atoms. 



The essence of our algebraic approach can be demon- 
strated using the following simple example. Let us con- 
sider a real symmetric 2x2 matrix, 



M{T) = 



El V 
V E2 



(1) 



which may represent the Hamiltonian of a two-state sys- 
tem. The notation implies that every matrix element 
is a function of the tunable parameter T, which could 
be an external field. The eigenvalues E of M{T), la- 
beled Ai and A2, may be found from the characteris- 
tic polynomial \M{T) - E\ = E^ + C\E + Co, where 
Co = E1E2 — V^ and Ci = —{Ei -I- -£2)- However, since 
the eigenvalues are themselves roots of this polynomial, 
i.e., {E — Xi){E — A2) = 0, the coefficients may equiva- 
lently be written in terms of the eigenvalues: Co = A1A2 
and Ci ~ — (Ai + A2). It should be noted that we did 
not explicitly calculate the eigenvalues to arrive at this 
conclusion. 

Since we are interested in curve crossings, we would 
like to know if the eigenvalues are degenerate, i.e., if 
Ai = A2, at some value of the parameter T. To this end, 
we introduce the discriminant A = (Ai — A2)^ of the 
characteristic polynomial of M{T) which can be written 
purely in terms of its coefficients : A = Cj — 4Co. If 
we now choose a simple parametrization for M{T), say 
£^1,2 = 2T and ^ = T, we find A = AT'^, a polynomial in 
T. The root of A, T = corresponds to the only crossing 
in the spectrum of M{T), as may be verified by explic- 
itly calculating the eigenvalues Ai^2 = T, 3T. The point 
of this example is to show how use of the discriminant 
maps the problem of finding spectral crossings to one of 
locating polynomial roots and bypasses the need to find 
eigenvalues. 

An extension of this simple example to larger matrices 
and more than one parameter leads to the study of mul- 
tivariate polynomials, a topic in algebraic geometry jj|. 
Assigning the task of a rigorous exposition to a separate 
publication [ig, we motivate the technique used in this 
Letter by analogy to the example stated above. We con- 



sider a matrix H{P), not trivially diagonal, of dimension 
n and depending polynomially on the set of parameters 
P = {Pi, ..Pn}- The characteristic polynomial of H{P) 
is of degree n in the eigenvalue E. The roots Ai.2,..n of 
this polynomial are used to define its discriminant 



D[H{p)] = n(^* 



A,f, 



(2) 



t<j 



which is an invariant of if (P) . It follows from Eq. Q that 
D[H{P)] does not change if the same quantity is added to 
every eigenvalue Xf, in fact it does not even change if the 
same quantity is added to each diagonal term of H{P), 
since this corresponds to a shift of the entire spectrum. 

Now as in the case of A above, D[H{P)] can also be 
calculated from the coefficients Ci of the characteristic 
polynomial without recourse to the roots. In the gen- 
eral case this can be done since D[H{P)] can be written 
in terms of the determinant of the {2n — 1) x (2n — 1) 
Sylvester matrix of the characteristic polynomial jj] : 



D[H{P)] 



(-1)- 



t^n 



On On — 1 ■•■ Oo 

Cn ■■■ Co 



Co 



nCn 



nC^ ... Ci 



2C2 Ci 



(3) 
Eq.® has two important implications. First, even if the 
roots cannot be calculated analytically, which is true in 
general for n > 5, D[H{P)] always can be. This means 
that the analytical mapping - supplied by the discrim- 
inant - between the level-crossing and polynomial root- 
finding problems is preserved for all n. Secondly, one 
can detect level crossings without solving for the spec- 
trum. This implies, as will be shown, the existence of 
powerful algebraic alternatives to numerically intensive 
spectral searches for crossings. 

As in the case of A the real roots of Eqn. 10) correspond 
to crossings in the spectrum of H{P); as well, the real 
parts of the complex roots of Eqn.ljSJ indicate avoided 
crossings jlTj . Further, in our simple example, A turned 
out to be quadratic in a single parameter T; in general 
D[H{P)] will be a high order polynomial in N variables. 
Typically, in order to locate the roots of this polynomial 
we will examine its coefficients. For example a simple 
method we will use is Descartes' Rule of Signs from ele- 



mentary algebra which equates the number of real roots 
of a polynomial to the number of sign changes in its coef- 
ficients, modulo 2. A more sophisticated concept we will 
use is that of a Sturm-Habicht sequence jj] of a polyno- 
mial, say in the variable Y. The difference in the number 
of sign variations in this sequence a.t Y ^ I and Y ^ r 
equals the number of real zeros of the polynomial in the 
interval [/, r] exactly. 

Let us now apply the methods outlined above to an 
atom in a uniform magnetic field. Such a system is rele- 
vant to the search for a Feshbach resonance, which typ- 
ically begins with a spin-polarized sample of atoms and 
requires a broad sweep of the magnetic field. Experimen- 
tally, if a crossing or a fine anti-crossing is encountered 
in the atomic spectrum while the magnetic field is tuned, 
some of the atoms can be transferred to a different spin 
state, introducing an impurity |l8l | . For theoretical book- 
keeping, it is useful to know if the low- and high-field 
states correlate adiabatically [l^- Such a correlation 
breaks down at a crossing. For these reasons, it is es- 
sential to be aware of crossing phenomena in the atomic 
spectrum. The magnetic fields typically used in Fesh- 
bach resonance searches are typically less than a Kilo- 
Gauss. Hence the parameter regime of our interest will 
heBFesh&[0,lkG]. 

To work with a concrete example we consider an ^^Na 
atom in its SS1/2 ground state split by the hyperfine in- 
teraction into two levels labeled by the angular momen- 
tum F = I + S where / = 3/2 and S = 1/2 are the 
nuclear and electronic spins respectively. In the presence 
of a uniform magnetic field B along the z-axis the atomic 
spectrum is described by the Breit-Rabi Hamiltonian |2flj 



Hbr - AI ■ S + B{aS, + bh 



(4) 



where A measures the strength of the hyperfine cou- 
pling and a and h may be easily related to the electronic 
and nuclear gyromagnetic ratios and the Bohr magneton. 
Throughout this Letter we will refer to the atomic eigen- 
states in the \F,AIp) basis where Mp is the component 
of F along the z-axis. In order to retain the polynomial 
form of the matrix elements, however, all representations 
will be made in the uncoupled basis \Ms,Mi), where 
Msj are the components of the respective spins along 
the z-axis. Such a representation in the ordered basis 

1=1 =3) |1 ^\ 1^ ^\ |1 ^\ 1^ 1\ |1 1\ 1^ 3) 
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J 



,3A - 2(a + 3b)E 



-3A + 2(a - 3b)E 
2V3A 



2V3A 
A - 2(a + b)£ 



-2(a - 
4A 



4A 

- 2[a 



A + 2(a + b)E 
2V3A 



2v^A 
-3A - 2(a - 3b)B 



3A + 2(a + 3b)BI 



To begin with a simplification we consider the fact that 
6 |<| a I in Eq.^ ^. Setting 6 = in Eq.© we cal- 
culate the discriminant of H'^ using built-in functions 
in Mathematica: 

J^i-^BRi — 65536 • ^ ' 

Considered as a polynomial in B, DlH^^] is even and 
exhibits no sign changes in its coefficients; by Descartes' 
rule it has no real roots for S 7^ |J|. Hence there are 
no level crossings in the spectrum. The real parts of the 
complex roots oi D[H'^] imply well documented, widely 
avoided crossings at S = and A/a (~ 316 G for ^^Na) 
[2l|. Physically, the latter corresponds to the crossover 
into the Paschen-Back regime. 

We now calculate the discriminant D[Hbr] for 6 7^ 0, 
the expression for which will be presented elsewhere on 
account of its length [la|. From the roots of D[Hbr] we 
infer the following changes introduced by the presence of 
b. The avoided crossings arc now located at S = and 
A/ {a — b) and have moved closer since b/a ~ —10"'^ for 
^^Na. More importantly there now appear crossings in 
the spectrum for B ^ 0. Evidently the coupling, b, of 
the nuclear spin to the magnetic field is the mechanism 
responsible for crossings in the Breit-Rabi spectrum. For 
B > there are 6 crossings clustered around 400 kG. In 
the regime Bpesh S[0, 1 kG] therefore the discriminant 
alerts us to the absence of any crossings as well as to the 
presence of two avoided crossings. 

To the best of our knowledge invariants of the Breit- 
Rabi Hamiltonian such as Eq.lO, which can be calcu- 
lated once the spins / and S are specified, have not been 
pointed out previously. They are complete catalogs of 
the parametric symmetries of Hbr- More generally, al- 
though we have picked an example from atomic physics, 
it should be clear that for any system represented by a 
parameter-dependent matrix, the discriminant is a pow- 
erful tool for investigating the location of crossings, their 
behavior as a function of the parameters as well as the 
physical mechanisms responsible for their occurrence. 

We now show how this method yields quite remarkable 
results in the more involved case of a diatomic molecule 
in a magnetic field. We consider a molecule made up of 
the previously considered ^^Na atom and an ^^Rb atom. 
The molecular Hamiltonian is 



H 



Mol 



K.E. + U{B) + V{R). 



(7) 



Here K.E. denotes the kinetic energy of the nuclei, U{B) 
the internal energies of the atoms in the magnetic field, 
and V{R) the molecular potential energy at the inter- 
nuclear distance R. Specifically, U{B) can be obtained 
from Eq.Ql, and 



V{R) = VsPs + VtPt, 



(8) 



where Vs.(^t) are the molecular (singlet, triplet) electro- 
static potentials, and Ps,t are projection operators on to 
the corresponding subspaces |19| . 



To work with a concrete example we prepare the Na-Rb 
atom pair in the state \2,2)Na\2,2)fii,, the "open" chan- 
nel for their collision. Since the component of the total 
angular momentum {Mptotai = 4) along the magnetic 
field is preserved by the interaction V(K)^ we need to 
consider in addition only the states which have the same 
Mptotai- |2,2}|3,2}, |2,1)|3,3), and |1,1)|3,3). Here we 
have retained the ordering of the atoms and dropped the 
subscripts. The four states mentioned above make up the 
electronic Hilbert space of the molecule. 

In the standard Born-Oppenheimer approximation 
(BOA) of molecular physics the K.E. in Eq.© is at 
first neglected and a set of adiabatic electronic potential 
curves obtained by diagonalizing the (in our case 4x4) 
interaction matrix 



H 



BO 



U{B) + V{R), 



(9) 



for a given B and for a large enough range of i?. In order 
to find the bound rovibrational levels one must then solve 
for the nuclear motion on these curves. The calculation 
has to be repeated for different values of B until a bound 
state equals the open channel in energy : this yields the 
simplest estimate of a Feshbach resonance location [l9| . 
It is imperative to know if there are crossings among 
the electronic potentials generated by Hboi because the 
BOA breaks down precisely at such points, and invali- 
dates the Feshbach resonance estimate. For this reason 
we investigate the discriminant D[Hbo]- But first we 
make a crucial transformation. We rewrite Eq.|(Sl) as 



V{R) ^ X{Ps - Pt) + {Vs + Vt)/2, 



(10) 



where X{R) — {Vs — Vt)/2, signifies the quantum me- 
chanical exchange energy responsible for chemical bond- 
ing and plays a critical role in the curve crossing prob- 
lem by allowing us to avoid the following difficulty. The 
Vs^T depend very sensitively on R. For ^^Na^^Rb for 
instance, the latest data requires 32(47) parameters to 
fit Vs[T) respectively [l^- These parameters are usually 
varied within their uncertainties, for instance to yield 
bounds on the Feshbach resonance predictions. It would 
require a numerical effort in (32-|-47-|-l=)80-dimensional 
parameter space, always limited by resolution, to inspect 
the spectrum for crossings for each Vs,Vt and B. We 
demonstrate below how the introduction of X renders the 
determination of curve crossings for Hbo completely in- 
sensitive to the complicated functional form of the Born- 
Oppenheimer potentials. 

Since D[Hbo] does not change if we subtract the same 
quantity from all the diagonal elements of Hbo i "we drop 
the second term in Ea. ((Tn|) . i.e. we calculate D[Hbo — 
{Vs + Vt)/2]{= D[Hbo])- We then substitute numbers 
for all the atomic constants in D[Hbo\ which gives us 
D[Hbo] = '^n=oPn{B)X" i.e., a polynomial of degree 6 
in X, whose coefficients depend on B. We now phrase 
our enquiry about curve crossings in the following way 



: is there any real value of the exchange energy {X £ 
[— oOj+cxd]) for which molecular curve crossings can be 
produced {D[Hbo] = 0) by the fields available in the 
laboratory [Bp^sh G[0,1 kG]) ? The answer can be found 
by examining the number of sign variations in the Sturni- 
Habicht sequence of D[Hbo] aX X ~ ±00 ^|. Only two 
real roots are found in the interval Bpesh G[0,1 kG], at 
S = OG and 502.2G, both for X = 0. Other than these 
we locate no crossings for any X and the BOA is valid 
everywhere else in the parameter regime of interest. 

We think it is quite remarkable, and points to the 
power of the algebraic approach, that in order to ex- 
tract information about crossings in the spectrum of a 
diatomic molecule in a magnetic field we did not avail of 
any knowledge about the Born-Oppenheimer potentials 
{ys,T) except that their difference is a real number! The 
curve crossings are entirely determined by atomic con- 
stants. Of course whether and at what R the crossings 
actually occur will depend on the specific form of the 
Vs,T- The algebraic method allows us to locate the com- 
plete set of magnetic fields that could cause the BOA to 
break down. If this set is not too restrictive, or is not in 
the neighborhood of a Fcshbach resonance, one need not 
refer to the Vs.t at all. 

Having demonstrated the application of curve-crossing 
methods to atomic and molecular spectra we now discuss 
possible extensions of our work. Our methods can readily 
be applied to other atoms and molecules. More generally 
we note that Eq.Q) is just a single example from an entire 
class of "spin-Hamiltonians" 



H, 



id 



y ^Bjjj. 



(11) 



Here the Ji are angular momentum operators and the 
indices i,j run over the spatial coordinates x,y,z. With 
the appropriate choices for the 'coupling' (C^) and 'field' 
(Si) parameters, effective Hamiltonians for many physi- 
cal systems can be set up and analyzed for level crossings. 
Specific examples of such systems are polar molecules in 
an electric field [22|, a spin-1 BEC with magnetic dipole 
interaction [23|, doubly-even nuclei described by the tri- 
axial rotor model |24J . massive spinors in an anisotropic 
universe |23, and so on. Of course our methods are not 
limited to Hamiltonian matrices. They can be applied 
to matrices that describe the stability of gyroscopes 3 , 
the polarization of light |2g , or the payoff in a game 3j , 
for instance. In sum, a large number of mathematical 
models contain parameters, and many of them involve 
matrices. We have shown how to extract curve crossing 
information from cases with effectively one and two tun- 
able parameters. Application to other systems as well as 
to a larger number of parameters are among the natural 
extensions of our work. 

In conclusion we have demonstrated a powerful and 
rigorous algebraic method for locating curve crossings 



in the spectra of physical systems. Along the way we 
have pointed out a new class of invariants of the Breit- 
Rabi equation of magnetic resonance and placed the 
Born-Oppcnheimer approximation for calculating Fcsh- 
bach resonances on a rigorous basis, using very little in- 
formation about molecular potentials. 
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